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Abstract 

With the help of transfer matrix method, the conditions for the ex- 
istence of the edge states in the semi-infinite armchair edged graphene 
is given. We discuss zero-energy and non-zero-energy edge states sep- 
arately, and show the nonexistence of the edge states in the model 
analytically and rigorously. 
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1 Introduction 

One of the interesting phenomena in solid state physics is the existence of 
edge states on the boundary, the properties of which are distinct from those 
of the bulk states, and they can play an important role in transport, there 
are examples showing that system is insulated in the bulk, while conduction 
can be allowed by edge states on the boundary. The most prominent ones are 
the quantum Hall effect (QHE) [H [21 131 S] and the quantum spin Hall effect 
(QSHE)[5l m [7], where the quantization of a Hall conductance is tightly 
associated with the edge states [21 EHIEIEIITIEIEIIIU]. 

The realization of graphene[TTl[T2l[T3 in laboratories has triggered great 
research interests in recent years. Prom the topological viewpoint, edge 
states can be induced in the system with different edge geometries [HI [15]; 
in experiments, with the help of scanning tunneling microscopy and scanning 
tunneling spectroscopy, the presence of structure-dependent edge states of 
graphite can be observed [Ml ITT]. There are some papers [181 [19] showing the 
similarity between the graphene model and the d-wave superconductor [20] 
with edges, where the existence of edge states depends on its edge shapes. 
For comparison, on a {110} surface of a dx2_y2 -wave superconductor, the 
origin of the Andreev bound states is the n -phase shift due to the uncon- 
ventional pairing symmetry, and graphene can be seem as an odd-parity 
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superconductor [19], at the zigzag/bearded edges, different signs of pairing 
potential lead to a vr-phase shift and zero-energy bound states, while at the 
armchair edge, it will not show the localized edge states. Hence the study 
on the zero mode of the edge state will have the fundamental meaning that 
can supply some information about the phase structure. 

The purpose of this paper is to study the edge states of the semi-infinite 
armchair edged graphene (AEG) analytically with the help of transfer matrix 
method [2H [221 ESI [Ml ES]- There are many papers showing the edge states 
in graphene models in different ways, while most of them focus on the zero 
mode of the zigzag graphene nanoribbon. And here we pay our attention on 
the semi-infinite armchair edged graphene, and study the zero-energy and 
non-zero-energy edge states. According to the discussions of the properties 
of the transfer matrix, we show an analytically proof of the non-existence of 
the edge states in the semi- infinite AEG. We have to mention that we just 
consider about nearest-neighbor-interaction in our whole discussions. 

2 EDGE STATES OF SEMI-INFINITE AEG 

The geometrical structure of graphene is shown in Fig{TJ It is infinite in y 
direction with periodic constant ^/3a, here a is the lattice constant between 
two n.n. A(B) atoms. We arrange all vertical chains in the order from the 
left to the right as {1, 2, 3, • • • }. Each vertical chain is an ID atomic chain, in 
which the periodic cell contains two kinds of atoms A and B. The position of 
the atoms are labeled by two indices {n,j), where n labels the order number 
of vertical chains from 1 to oo, with n = 1 representing the left edge. 

Due to its infinity in y direction, ky is a good quantum number, so 

that we can take the Fourier transformation on wave functions ^^('^■')(^?/) 

for A{B) atoms: = ^w{'''kyyj^A{B))^n.A{B){f^y)- The Hamiltonian 

now in {n, ky} representation can be expressed by a set of Fermion operators 
{^n,A(B){ky)},n = 1, 2, . . . , oo: 

ky,n^l 

+ t e-^'y'^/'^'^ljky)^r.-^Mky) 
+ t Y e-^'y^/'^'^ljky)^n+iMky)+ h.c. (1) 
Prom the dynamical equations of {$„^^(^)(A;y)}: E'<l>„^^(^)(/cy) = [<I>„^^(5)(/cy), ff]. 
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Figure 1: Schematic illustration of the lattice structure of semi- infinite arm- 
chair edged graphene. The black circles represent the type B atom, while 
empty circle represents type A atom; It is infinite in y direction and the 
edge locates at the left perpendicular chain labeled by number 1. 



we get following recursive equations for each fixed wave vector ky: 



(2) 



where we use ^n,A{B) instead of ^n,A{B)if^y)- The equations of the wave 
functions for the boundary sites (n = 1) are 



E<^i^A 



(3) 



We can get the bulk energy for the extended states, according to Eq. (I2|), 
take kx as good quantum number, the bulk energy satisfies: 



E'^ = 1 + 4cos^ {kxa/2) + 4cos {V3kya/2) cos {kxa/2) 



(4) 



which is the energy dispersion relation of the 2D graphene system. In the 
following, we will focus on the edge states, and discuss the zero-energy {E = 
0) and non-zero-energy (E ^ 0) cases, separately. To examine whether the 
edge states can exist in the semi- infinite graphene model or not. 
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(1) E = 0: the zero-energy spectrum. The recursions Eq. ([2]) are reduced 
into two decoupled equations for sublattice A and B: 



a = A,B, r]A = 1, f?B = -1, n = 2, 3, 



(5) 



Above equations can be rewritten as follows transfer matrix equations: 



l.Q 5 



n = 1,2,3,- 



where 



-tr), 



,%/3kya/2 



(6) 
(7) 



_g-iriaV3kya/2 

1 



(8) 



detT„ 



A 



(a) 



1. It is easy to get two eigenvalues of Ta{a = A,B): a|^^ = 
^-zriakyaV3/2 ^ |Aa| = 1 and two eigenvalues 
(A^)*. It results |A^| = 1, due to A" -A^ = 1. 



^Aq, ± liy^4: - Al where Aq 



_i_ must be complex and A^ 



Thus, we can conclude that the states are always extended in the semi- 
infinite AEG, no zero-energy edge state can exist. 

(2) E ^ 0: according to Eq. ([2]), the equations of two sublattices are 
coupled with each other. We can also write down a transfer matrix for A 
or B type sublattices, and it becomes 4x4. For more concise description, 
we define a 4 x 1 vector its transpose (^^^)* is a 1 x 4 vector as 



follows: 



n—l,a 



n—2,a 



^n-3,a ) • Then, the Eq. ([2]) can be expressed as 



^ a 



( G \ 

1 

1 





(n+l) 
a , 

G 


1 



n ^ 3 
-1 \ 





/ 



det To = 1 



(9) 



where G = -26, 6 = cos [V3kya/2), A = {E'^-3t^)/t^ = E'^-3 (set t = 1). 
We further introduce one fictitious lattice line at left of edge chain and set 
them all to be zero: ^ = 0. It will not change anything even introducing 
the hopping coupling with 1st line. Now Eq. ^ and Eq. ([3]) can be written 
as follows: 



^ a 
$(3) 



rpn- 

^0 



G*(,2) + ^ 



(1) 

a • 



(10) 

(11) 
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the eigenvalues of transfer matrix Tq can be written in the form of {Aj}: 

l) , A2 = ■ 



Ai 
A3 



(^-02 - ^/aj-A^ , A4 = i (^-02 + \/a2 - 4 



(12) 



where 

ai = 6* + ^, 02 = 6'-\/m, m = A + 2 + 6l2. (13) 

Now, define a transformation {/, where exists, to rewrite the transfer 
matrix and its relation Eq. Q with D = U^^TqU = diag{Xi, X2, X3, 
In order to get the physically meaningful states, we have to ensure that all 
eigenvalues {|Ai|} must be finite for any finite energy E. 

When m 7^ 0, of 7^ 4, and 7^ 4, we can write down U as follows 



U 



( Af A| A| A^ ^ 
Ai A2 A3 A4 
Ai A2 A3 A4 

Villi/ 



(14) 



with det [/ 7^ 0. Its inverse matrix U ^ can be found and its elements of 
U^^ are described as follows 





= Cij/ui{ai - 02), 








Ui = 


af - 4 + ai^a'f - 


4 = 


-2AiVaf 


-4 


< U2 = 


A-aj + aiy/a'f- 


4 = 


-1\2^d\ 


- 4 


U3 = 


a^- 4. + a2^d\ - 


4 = 


-IX^^/ai 


- 4 


^Ui = 


02 — 4 — a2\/a| — 


4 = 


IXi^J CL^ — 


4. 



(15) 



where 



C41 
C12 
C22 

Ci3 

C23 
C33 
C43 
C32 
C42 

Ci4 

C34 



31 



C21 — — Cii — 2, 
- 4 - ai - 2a2, 



\/af - 4 + ai + 2a2 
-2 - 0102 + a2-\/ai - 4 
2 + aia2 + a2 \/ai - 4 
2 + aia2 



aiVc 
2 + aia2 + ai^/c 



4, 



(16) 



2ai + 02 - \/a| - 4, 



2ai + 02 + \/a2 - 4, 
-ai + 
012 - 



- 4, C24 = ai + ^of 
4, C44 = 02 + vof- 
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The relation JJ-^TqU = D was carefully confirmed by U and U 
tedious but straightforward algebraic calculations, we obtain 



^(n+3) ^ ;7/)"^-l^(3) ^ jj 



W2X2 
\ W4X2 ) 



,n > 1. 



After 



(17) 



Where 



{U-2' + GU-,') ^2,a + (^-3' + (A + 1) U-,') 

i = 1,2,3,4. 



(18) 



Since all elements of matrix Uij 



are non-zero and finite, thus the criteria 
whether the state $n,a is extended or edge state can be determined by those 
{H^jA"} regardless the front matrix U in Eq. ()17p . 

From the definition of {Aj}, we can confirm Ai • A2 = 1 and A3 • A4 = 1. 
thus IA2I = |Ai|~^ and IA3I = |A4|~^. Define {Ai,A2} ({A3,A4}) as a pair 



partner. A 



2(1) 



is the partner of Ai(2)i while A4(3) 



the partner of A3(4) . Then 



Xj and Ai 



A2. It 



we use A^ to denote the partner of Aj. For example, A2 
is similar to define the partner of Wi [see Eq. pzp and Eq. (|18|) ]: the partner 
of PFi(2) is VF2(i) = Wj(2), and W3(4) the ^4(3) = W^jy If |Ai| = 1, we do 
have = 1, which corresponds to the extended states. As we know, the 
edge states correspond to some |Ai| < 1, due to the properties of Aj and its 
partner Xj, in order to get the physically meaningful non-zero-energy edge 
states, the amplitude of which must be finite at n — )• 00 in Eq. ()17p . so that 
we can get the necessary condition for the edge states Wi = and Wj 7^ 
when |Ai| > 1 and \Xj\ < 1. 

Without loss of generality, at first we assume |Ai| > 1 at some {E, ky}. 
In the case, we have IA2I < 1- The necessary condition for the existence 
of edge state is Wi = and W2 = rj ^ 0. From Eq. (fT8]) . the necessary 
condition turns to 



Where 



All^2,a + ^12^1,a = 0, 
A21^2,a + A22^l,a=r]{^ 0). 



(19) 



Au = U^2 + GU^i^ Ai2 = t/iV + (A + l)U^i 



T-1 



-1 



r-1 



A2i = U^i + GU:^,\ ^22 = ^/2-3' + (A + l)f/2l' 
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In terms of relations ([TS]) , (fT6]) and (fT3|) , we get 



^11 



Ci2 + GCii 1 



^12 



Mi(ai-a2) 2i/af - 

Cl3 + (A + l)Cii _ 02 



^21 



ni(ai-a2) 2Y^af^\/m' 
C22 + GC21 —1 



^22 



M2(ai - 02) 2v^af - A^/m 
_ C23 + (A + l)C72i _ -02 



(20) 



^2(01- 02) 2Y^af - A^/m' 



where the relations Eq. (fT2]) . Eq. (fT3]) . Eq. (fTSj) . and Eq. (fT6|) have been 
used. According to m 7^ and af 7^ 4 and a| 7^ 4, Au and y42i both are 
finite and not zero. If 02 = 0, we have A12 = A22 = 0. The 1st Eq. ()19p 
must lead to ^2.a = because of An 7^ 0, meanwhile the 2nd Eq. ()19p 
showing ^22'l'i,a = V ^ could not be satisfied. Thus 02 = must result 
W2 = 0. Assume all elements {^jj} are not zero (02 7^ 0) in our following 
discussion. To ensure that both equations in Eq. (fT9j) should be satisfied 
with r/ 7^ 0, we can write down a non-homogeneous linear equations with 
vector {^2,a, ^i,aY, therefore Eq. (fT9]) can be written as 



where ^ is a matrix with the elements {Aij : i,j = 1,2}, and suppose 
b = (0, tjY- The existence of the solution above needs rank(A) = rank(^, b). 
After simple but tedious algebraic calculation, we can find that rank(^) 7^ 
rank(A, 6), when r] ^ 0. So we can conclude that there is no non-zero-energy 
edge state here. Besides, we can get similar relations with other Aj > 1 and 
Xj < 1. The details are neglected here, all the discussion are rigorous. Thus 
we rigorously proved that no any edge state exists. 

When m = 0, that leads to ai = 02 = 9, {A,} are still the eigenvalues of 
transfer matrix Tq. Because the values of 1st and 3rd column of U are the 
same, and those of 2nd and 4t/i the same as well, it is easy to get det U = 0, 
so doesn't exist. And it is easy to confirm that |Ai| = IA2I = 1 and 
I A3 1 = IA4I = 1, which corresponds to the oscillation of the wave function, 
without decay. Meanwhile, under this condition, we can get E'^ = \ — (P'. 
Taking this value back to Eq. we can get the relation 2 cos [kxa/2) = 
cos (\/3fcyo/2). Energies show in Eq. @ are complete with both real wave 
vectors, so it results that there are just bulk states when m = 0. 

When af = 4 or = 4, does not exist, because the values of 1st and 
2nd column of U are the same. Without loss of generality, we can choose 
af = 4 in our following discussion, and it is easy to know m is real, so we 
have oi = 2. Then we can easily get |Ai| = IA2I = 1, while at some ky, 




(21) 
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I -^3 1 = I -^4 1 = 1 corresponds to the extended states, and others IA3I < 1 and 
IA4I > 1 just A3-A4 = 1. Although doesn't exist, we can suppose to write 
down it with the definition of right inverse matrix at ai — )• 2, so we can find 
at this hmitation, the 1st and 2nd row of are closely the same in values, 
while both of them correspond to the extended states. And besides, we can 
consider the 3rd and 4:th row which corresponds to the possible decay wave 
with a certain ky . Then it is easy to get similar relation as Eq. p8|) . In order 
to get the physically meaningful edge states, we can get similar relations as 
Eq. p9|) . Considering about a similar inhomogeneous equations, we can 
conclude that there is no non-zero-energy edge state under this condition. 
Besides, we can give the similar discussion with 02 = —2. 

In this paper, we give an analytically proof of the non-existence of the 
edge states in the semi-infinite AEG. By this method, we show the bulk 
energy spectrum and the condition for the existence of the physically mean- 
ingful edge states, according to the rigorous discussion of the transfer matrix 
and its eigenstates, we can get the necessary and sufficient condition for the 
existence of zero-energy and non-zero-energy edge states, and finally we find 
the contradictory condition to show the non-existence of them. 
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